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We present a 2.5<r detection of the Integrated Sachs- Wolfe (ISW) effect and discuss the constraints 
it places on cosmological parameters. We cross-correlate microwave temperature maps from the 
WMAP satellite with a 4000 deg 2 luminous red galaxy (LRG) overdensity map measured by the 
Sloan Digital Sky Survey. These galaxies have accurate photometric redshifts (Az ~ 0.03) and an 
approximately volume limited redshift distribution from z ~ 0.2 to z ~ 0.6 well suited to detecting 
the ISW effect. Accurate photometric redshifts allow us to perform a reliable auto-correlation 
analysis of the LRGs, eliminating the uncertainty in the galaxy bias, and combined with cross 
correlation signal, constrains cosmological parameters - in particular, the matter density. We use a 
minimum variance power spectrum estimator that optimally weights the data according to expected 
theoretical templates. We find a 2.5u signal in the Ka, Q, V, and W WMAP bands, after combining 
the information from multipoles 2 < I < 400. This is consistent with the expected amplitude 
of the ISW effect, but requires a lower matter density than is usually assumed: the amplitude, 
parametrized by the galaxy bias assuming Qm = 0.3, Qa = 0.7 and erg = 0.9, is b g — 4.05 ± 1.54 
for V band, with similar results for the other bands. This should be compared to b g = 1.82 ± 0.02 
from the auto-correlation analysis. These data provide only a weak confirmation (2.5<r) of dark 
energy, but provide a significant upper limit: Qa = 0.80lo'gg (lcr)^ ) ' ( j'g(2cr), assuming a cosmology 
with Qm + Qa = 1, Qb = 0.05, and erg = 0.9, and w = —1. The weak cross-correlation signal rules 
out low matter density/high dark energy density universes and, in combination with other data, 
strongly constrains models with w < —1.3. We provide a simple prescription to incorporate these 
constraints into cosmological parameter estimation methods for (Qm, erg, w). We find no evidence 
for a systematic contamination of ISW signal, either from Galactic or extragalactic sources, but we 
do detect some large statistical fluctuations on smaller scales that could affect analyses without the 
template weighting. 

PACS numbers: 



I. INTRODUCTION 

The cosmic microwave background (CMB) observed by 
the WMAP satellite pj has been a font of information 
for cosmology. The positions of the acoustic oscillations 
in the angular power spectrum of temperature fluctua- 
tions imprinted at the last scattering surface measures 
the curvature of the universe to unprecedented accuracy. 
In addition, a careful modelling of the power spectrum 
combined with other cosmological probes have allowed 
extremely precise measurements of the parameters of the 
ACDM model. Subsequent measurements of the power 
spectrum by WMAP , as well as future experiments such 
as the PLANCK satellite will further constrain the space 
of cosmological models, as well as test our understanding 
of the physics underlying recombination. 

CMB temperature fluctuations are not just sourced 
by density fluctuations at the last scattering surface; 
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anisotropics also arise due to the interaction of photons 
with hot electrons in galaxies and clusters (the Sunyaev- 
Zeldovich [SZ, and kinetic-SZ effects), as well as 
with the gravitational potentials along their propaga- 
tion path (the Integrated Sachs- Wolfe [ISW,|j] effect and 
gravitational lensing) . 

The ISW effect results from the red- (or blue-)shifting 
of CMB photons as they propagate through gravitational 
potential wells. If these potentials did not evolve, then 
the blueshift gained falling into a potential well would ex- 
actly cancel the redshift emerging from the well; evolving 
potentials spoil this cancellation. This is significant only 
on the largest scales, both because the power in the fluc- 
tuations of the potential is largest on large scales, and 
because integrating along the line of sight cancels out 
the temperature fluctuations (a photon is as likely to be 
redshifted as blueshifted) . 

Although the above effects are too small to unam- 
bigously detect in current CMB data (for eg. |5|), 
one can attempt to isolate them by cross-correlating 
the temperature maps with suitable tracer populations 
HHBSHmillEi- This P a P er is a continuation 
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of these efforts, attempting to detect the ISW effect, but 
also goes one step further in making the ISW signal use- 
ful for cosmological parameter analysis. We combine the 
cross-correlation analysis with a galaxy auto-correlation 
measurement to remove the uncertainty due to galaxy 
bias. The auto-correlation analysis is only reliable when 
the galaxy population has well determined redshifts: we 
argue below that this is indeed the case for the sam- 
ple used here. This paper does not focus on the auto- 
correlation analysis, so we perform a restricted analysis 
where we assume the shape of power spectrum and only 
constrain the amplitude. This is sufficient for the cur- 
rent purposes, as the errors are dominated by the sta- 
tistical uncertainties in the ISW effect. A detailed auto- 
correlation analysis will be presented in a future paper. 

A galaxy catalog used for detecting the ISW effect 
must satisfy several criteria. Since the ISW signal is only 
detectable on the largest scales, the catalog must cover 
as large an area as possible. In addition, a large num- 
ber of galaxies is necessary to remove noise from random 
Poisson fluctuations. Finally, to make theoretical inter- 
pretation possible, these galaxies must be drawn from 
a uniform population with a well characterized redshift 
distribution. 

Photometrically selected luminous red galaxies (LRGs) 
H3 from the Sloan Digital Sky Survey (SDSS) [3] are 
a promising candidate. These galaxies are amongst the 
most luminous in the Universe and so probe cosmolog- 
ically interesting volumes, including the transition from 
matter to dark energy domination that sources the ISW 
signal. These galaxies are old stellar populations with 
very uniform spectral energy distributions that make uni- 
form selection and accurate photometric redshifts possi- 
ble. Finally, the deep and wide-angle imaging of the sky 
by the SDSS allows these galaxies to be selected over a 
large area with high densities. 

Our goal is to extract the maximal amount of informa- 
tion on the ISW effect from the data. To do so, we employ 
quadratic estimator methods, which combine the data in 
an optimal way in the presence of noise and incomplete 
sky coverage In addition, we weight the infor- 

mation on different scales according to the expected scale 
dependence of the signal: the ISW effect is expected to 
dominate on large scales and be absent on small scales, 
so we weight these different scales accordingly. We ap- 
ply these techniques to the latest SDSS data reductions, 
ensuring the largest available SDSS sky coverage to date. 

Detecting an ISW effect would provide evidence of ei- 
ther a cosmological constant, quintessence or curvature 
independent of supernovae observations of the luminos- 
ity distance 0. We investigate the cosmological impli- 
cations of an ISW detection within the assumption of a 
flat universe, and in particular, on its ability to constrain 
the matter density. The rapid increase of the ISW signal 
with decreasing matter density allows us to place strong 
lower bounds on the matter density and suggests an al- 
ternative probe of the properties of the dark energy. 

The paper is organized as follows : Sec. [H] intro- 



duces the ISW effect, developing the formalism to com- 
pute the predicted ISW signal. Sec. 11111 describes the 
WMAP and SDSS data used in this paper. Sec.HVland 
Sec. then discuss the cross correlation of these data, 
and the results obtained. Systematic effects are discussed 
in Sec. lVll while cosmological implications are considered 
in Sec. IVIll We conclude in Sec. IVIlH Unless otherwise 
specified, we assume a ACDM cosmology with Qm = 0.3, 
Q A = 0.7, H = lOO/i km/s/Mpc and <j 8 = 0.9. Where 
necessary, we use h — 0.7. 

II. THEORY 

We briefly review the ISW effect, and its cross- 
correlation with the galaxy density (see also eg. [|| 0, 
113, El) We consider linear functions of the density field 
projected onto the sky, specializing to the cases of the 
ISW temperature perturbations and the galaxy overden- 
sity field. Decomposing these projected fields into spher- 
ical harmonics allows us to compute the expected auto- 
and cross-power spectra. For simplicity, we restrict our 
discussion to flat universes (f2 m + tt\ = 1 ). 

A. Projections onto the sky 

We start with the 3D matter overdensity field, £30 (y), 
assumed to be an isotropic random variable. This allows 
us to define the 3D matter power spectrum, P(k), 

<<Mk)Wk')> = (2^) 3 <5(k - k')P(fc) • (I) 
We project this density field onto the sky, 

p(n) = J dyf(y)£[S 3D (y,yh)} , (2) 

where f(y) weights the density field as a function of co- 
moving distance, y, and £ is a linear operator (indepen- 
dent of y) operating on the 3D density field. Simulta- 
neously expanding the projected field in spherical har- 
monics, Y; m (n), and Fourier transforming the 3D density 
field, we obtain 

Plm=i l j J dyfiyMky^MClSMvM > 

(3) 

where we use the orthogonality of the Yi m 's and the ex- 
pansion of a Fourier wave in spherical coordinates, 

e-^y = 4* £ *^(fcy)^„(k)y ; * m (n) . (4) 

I in 

1. Galaxy and Temperature Fluctuations 

Given the above equation, we specialize to the galaxy 
overdensity and the ISW temperature fluctuations. We 
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assume that the galaxy overdensity, S 9t3 D, is related to 
the matter overdensity by a linear bias, b g , 

5g,3D = £g[<>3D] = b g S 3D , (5) 

which implies that (Eq. |3j) , 

d 3 k 



5g,i m = i l b g / / d v f(y)ji(ky)Yi* m (k)5 3D {y,k) 



(6) 

Although linear bias must break down at some scale, it 
appears to be well motivated both theoretically!^ and 
observationally on the large scales that the ISW ef- 
fect is sensitive to. The redshift weighting f(y) is simply 
given by the galaxy selection function, </>(y), appropri- 
ately normalized, 



m = 



y 2 0(y) 

J dyy 2 4>{y) 



Working on large scales allows us to use linear perturba- 
tion theory where the growth of density fluctuations with 
time is separable, 



S 3D (y,k) = D(y)S 3D (k) , 



(8) 



where D(y) is the growth factor. Substituting this above 
allows us to write Eq. m the following useful form, 

/d 3 k 
— 5 3D (m*MW g (k) , (9) 

where the window function W g {k) is given by 

W g (k)= J dy f(y)D(y)j l (ky) . (10) 

The temperature fluctuations due to the ISW effect 
are given by the line of sight integral of the change in the 
gravitational potential to the last scattering surface, 



/AT 



ISW 



ii'i 



= -2 / dy$(y,yn) 



(11) 



where $ is the gravitational potential, and the overdot 
is the derivative with respect to conformal distance (or 
equivalently, conformal lookback time) at constant y 3 D- 
The gravitational potential can be related to the 3D den- 
sity fluctuations by Poisson's equation, 



v 2 $ = -met 



'ZD 



2~ o"m— - , 
where a is the scale factor, implying that, 



(12) 



d 



dy 



C ISW [S 3D ]^-3H^ M i-[y- 2 —) ■ (13) 



a 



Taking the Fourier transform and substituting into Eq.|21 
we find that, 



/AT 



ISW,lm 



/ d 3 k 3n M H$ , 

2vr 2 k 2 ' ay 



^{ky)Y^)l. (5Mv ^ 



dy V a{y) 



(14) 



Again restricting to linear theory allows us to write the 
above equation in the same form as Eq. [5J 



/AT 



ISWlm 



d 3 k 

(15) 



where the ISW window function is given by 

30 H 2 f Vo / D\ 

W ISW (k) = ^p- J dvMky)dJ^y (16) 

It is straightforward to compute the relevant power 
spectra. The galaxy power spectrum, Cf 9 = (6g t i m S* lm ) , 
is given by, 



Cf 4 ,T / dk^\Wg(k)\ 



(17) 



where we have used the definition of the three dimen- 
sional power spectrum (Eq. ^) and A 2 (k) is the variance 
per logarithmic wavenumber, 



1 

(2~7) 3 



A 2 (fc) = — ^-r47rfc 3 P(fc) 



(18) 



The cross correlation between the ISW temperature fluc- 
tuations and galaxy overdensities is similarly given by, 



g-ISW 
I 



A 2 (k) 

An I dk -^lw g {k)W ISW {k) . (19) 



B. Predictions 

The above formalism allows us to compute the ex- 
pected ISW power spectrum given the redshift distribu- 
tion of a tracer population and a particular cosmology. 
Fig.^shows the ISW power spectrum assuming the LRG 
redshift distribution introduced below (Sec. 1111 B|) and a 
flat cosmology. The growth factor is given by 



D(a) 



5 H{a) 

n"M,CT 



da' 



H Jo {a'H{a')/Hof ' 



(20) 



where a = 1/(1 + z) is the scale factor, and H (a) is deter- 
mined by the Fricdmann equation. We note that (within 
the context of flat cosmologies), D/a is independent of 
redshift if CIm = 1, leading to the lore that the detec- 
tion of the ISW effect is evidence for some form of dark 
energy, although we note that universes with curvature 
and no dark energy /cosmological constant also lead to 
an ISW effect. 

Finally, a computational note on performing the spher- 
ical bessel integrals of the previous section : we follow 
[24| and recast these as logarithmically discretized Hankcl 
transforms. In this form, the integrals can be efficiently 
performed via FFT convolutions using the FFTLOG al- 
gorithm [25| as implemented in |24j. 
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FIG. 1: Predictions for the ISW signal given the redshift 
distribution of LRGs described in Sec. 1111 51 b g = 1, as = 0.9, 
and a flat universe. The different curves show the effect of 
changing the matter density. In particular, observe that the 
effect becomes stronger as the matter density decreases. 
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FIG. 2: The LRG redshift distribution. The histogram shows 
the photometric redshift distribution, the curve is the true 
redshift distribution estimated by regularized deconvolution 
of the photo- z errors. The dotted lines show the photometric 
redshift cuts imposed at z = 0.2 and 0.6. 



III. DATA 
A. CMB temperature from WMAP 

The WMAP mission pj is designed to produce all-sky 
maps of the CMB at multipoles up to I ^several hundred. 
This analysis uses the first public release of WMAP data, 
consisting of the first year of observations from the Sun- 
Earth L2 Lagrange point. WMAP carries ten differencing 
assemblies (DAs), each of which measures the difference 
in intensity of the CMB at two points on the sky; a CMB 
map is constructed from these temperature differences as 
the satellite rotates. WMAP observes in the K (1 DA), 
Ka (1 DA), Q (2 DAs), V (2 DAs), and W (4 DAs) bands 
corresponding to central frequencies of 23, 33, 41, 61, and 
94 GHz, respectively. The WMAP team has pixelized the 
data from each DA in the HEALPix |6lJ pixelization sys- 
tem at resolution 9 [HS^]. This system has 3,145,728 
pixels, each 47.2 sq. arcmin in area. These maps are not 
beam-deconvolved; this, combined with the WMAP scan 
strategy, results in nearly uncorrelated Gaussian uncer- 
tainties on the temperature in each pixel. 

We limit this analysis to the Ka through W bands. The 
K band is heavily contaminated by Galactic emission, 
increasing the number of possible systematic effects and 
making any error analysis unreliable. We apply the KpO 
mask specified by the WMAP team to mask out regions 
where Galactic foregrounds dominate. In addition, we 
reject pixels in the WMAP point source mask, including 
a 0.6 degree exclusion radius around each source. Finally, 
since the SDSS area used only covers about a tenth of 



the sky, we eliminate all pixels greater than 10 degrees 
from the SDSS mask, allowing us to speed up our cross- 
correlation analysis. This mask, denoted Kp0nS10\ps , 
admits 756,078 HEALPix pixels (9915 sq. deg). The 
variance in each of these pixels is computed directly from 
the number of observations of each pixel, as specified by 
the WMAP team. 

We chose not to use either the WMAP "Internal Linear 
Combination" (LLC) map j2||, or the foreground cleaned 
map of |28|. to avoid a number of practical difficulties. 
These maps lose the frequency dependence of the original 
maps (useful to identify contaminating signals) and have 
complicated pixel-pixel noise correlations (increasing the 
complexity of the error analysis). 

B. SDSS Luminous Red Galaxies 

The Sloan Digital Sky Survey is an ongoing effort 
to image approximately ir stcradians of the sky, and ob- 
tain spectra of approximately one million of the detected 
objects [2^, The imaging is carried out by drift- 

scanning the sky in photometric conditions |3lj |. in five 
bands (uari z) l32t |3jj using a specially designed wide- 
field camera |3J|. Using these data, objects are targeted 
for spectroscopy and are observed with a 640-fiber 
spectrograph on the same telescope. All of these data 
are processed by completely automated pipelines that de- 
tect and measure photometric properties of objects, and 
astrometrically calibrate the data [36l l37j . The SDSS 
is nearing completion, and has had four major data re- 
leases [33, HH, [H^] ■ This paper uses all data observed 
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th roug h Fall 2003, processed and calibrated as described 
in HJ. 

The usefulness of LRGs as a cosmolo gica l pr obe has 
been appreciated by a number of authors [30L |42| . These 
are typically the most luminous galaxies in the universe, 
and therefore probe cosmologically interesting volumes. 
In addition, these galaxies are generically old stellar sys- 
tems with uniform spectral energy distributions (SEDs) 
characterized principally by a strong discontinuity at 
4000 A. This combination of a uniform SED and a strong 
4000 A break make LRGs an ideal candidate for pho- 
tometric redshift algorithms, with redshift accuracies of 
o z ~ 0.03 Q. We briefly outline the construction of the 
photometric LRG sample used in this paper below; a de- 
tailed discussion of the selection criteria and properties 
of the sample is in [3 . 

Our selection criteria are based on the spectroscopic 
selection of LRGs described in [j^j , extended to lower ap- 
parent luminosities. We select LRGs by choosing galaxies 
that both have colors consistent with an old stellar pop- 
ulation, as well as absolute luminosities greater than a 
chosen threshold. The first criterion is simple to imple- 
ment since the uniform SEDs of LRGs imply that they lie 
on an extremely tight locus in the space of galaxy colors; 
we simply select all galaxies that lie close to that locus. 
More specifically, we can define three (not independent) 
colors that describe this locus, 

c± = (r - i) - 0.25(5 - r) - 0.18 , 

d± = (r -<)- 0.125(g -r) , 

C |j = 0.7(5 - r) + 1.2(r - i - 0.18) , (21) 

where g, r, and i are the SDSS model magnitudes [3^ 
in these bands respectively. We now make the following 
color selections, 

Cut I: | CJ J<0.2; 
Cut II : d± > 0.55, g - r > 1.4. (22) 

Making two cuts (Cut I and Cut II) is convenient since 
the LRG color locus changes direction sharply as the 4000 
A break redshifts from the g to the r band; this division 
divides the sample into low redshift (Cut I, z < 0.4) and 
high redshift (Cut II, z > 0.4) samples. 

In order to implement the absolute magnitude cut, 
we follow [3(j and impose a cut in the galaxy color- 
magnitude space. The specific cuts we use are 

Cl 

Cut I : r Petro < 13.6 + — , r Petro < 19.7, 
Cut II: i < 18.3 + 2d ± , i < 20, (23) 

where rp etro is the SDSS r band Petrosian magnitude 
|S8j . Finally, we reject all objects that resemble the point- 
spread function, or if they have colors inconsistent with 
normal galaxies; these cuts attempt to remove interloping 
stars. 

Applying these selection criteria to the ~ 5500 degress 
of photometric SDSS imaging in the Galactic North 
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FIG. 3: The LRG angular distribution in Galactic coordi- 
nates. The gaps in the distribution are due to the stellar 
mask, nonphotometric data, and Galactic extinction. 



yields a catalog of approximately 900,000 galaxies. Ap- 
plying the simple template fitting photometric redshift 
algorithm of [lj| > we restrict this catalog to galaxies with 
0.2 < z photo < 0.6, leaving us with ~ 650,000 galaxies. 
We use the regularized inversion method of [14J as well 
as the photometric redshift error distribution presented 
there, to estimate the true redshift distribution of the 
sample. The results, comparing the photometric and true 
redshift distributions are shown in Fig. The LRGs are 
approximately volume limited from z ~ 0.2 to z ~ 0.55. 
Comparing this to the S /N estimates of 0] > we see that 
the LRG redshift distribution is well suited to detecting 
the ISW effect. Therefore, despite the availability of ac- 
curate photometric redshifts, we do not weight the LRGs 
in redshift any differently than they already are in the 
sample we use here. 



This catalog is pixelized as a number overdensity, 
S g = Sn/n, onto a HEALPix pixelization of the sphere, 
with 3,145,728 pixels. Since the photometric catalogs are 
incomplete around bright stars, we mask regions around 
the stars in the Tycho astrometric catalog |45|. We also 
exclude data from the three southern SDSS stripes due 
to difficulties in photometrically calibrating them rela- 
tive to the data in the Northern Galactic Cap. The 
final catalog covers an solid angle of 3,893 square de- 
grees (296,872 HEALPix resolution 9 pixels) and contains 
503,944 galaxies at a mean density of 1.70 galaxies per 
pixel. The sky coverage is shown in Fig. [3] 
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IV. ESTIMATING THE CROSS CORRELATION 

We start by organizing the temperature fluctuations 
and the galaxy overdensities into a single data vector, 



x = (x BjT ,x ff ) , 



(24) 



where xg^ is a vector with the measured CMB tem- 
perature (with the monopole and dipole subtracted) in 
band B at every HEALPix pixel; analogously, x 5 is 
the LRG number overdensity. We suppress the band 
subscript for simplicity, with the implicit understanding 
that we always refer to the cross correlation of a single 
WMAP band with the LRG overdensity. The covariance 
matrix of x is, 



where C dlag is given by, 



C 9l 



C 9T 



&diaq — 



C TT + N TT 



C99 + N99 



(25) 



(26) 



where N xx is the pixel noise matrix. The submatrices 
C TT , C 39 and C gT are defined by 



a 



ab 



l in 



(27) 



where nf is the position (on the sky) of the i th point of the 
vector x a . The temperature-temperature, galaxy-galaxy 
and galaxy-temperature angular power spectra are de- 
noted by Cj T , Cf 9 and Cf T respectively. In our analy- 
sis below, we use the best fit ACDM prediction for the 
temperature angular power spectrum, multiplied by the 
appropriate WMAP instrumental response at each mul- 
tipole. The galaxy power spectrum is estimated using a 
pseudo-C; estimator ^3| , and fit by the non- linear power 
spectrum of |46|. multiplied by a constant linear bias. We 
project out the monopole and dipole of both these power 
spectra by setting the power in the I = 0, 1 modes to 
a value (10 -2 ) much greater than the true power spec- 
trum. In addition, the measured galaxy-galaxy power 
spectrum has excess power on large scales (I < 20) due 
to calibration errors; we boost the prior power spectrum 
to account for this. 

We parametrize Cf T as a sum of bandpowers, P,,;, with 
amplitudes c, to be estimated, 



c 



gT 



(28) 



We consider two bandpowers in this paper. The first are 
"flat" bandpowers given by 



Pi,l — -^(0 s li.rnin ^1^1% 

— ; otherwise , 



(29) 



where B{1) is the WMAP instrumental response |47j. 
This parametrizes the power spectrum as a sum of step 



functions and is useful when the shape of the power spec- 
trum is unknown. However, as discussed in Sec. [nj the 
shape of the ISW component of the galaxy-temperature 
correlation is well determined by the cosmology and the 
redshift distribution of the galaxies with only its ampli- 
tude (equivalently, the galaxy bias) unknown. We define 
the "template" bandpowers by replacing the flat band- 
powers for I < l max with the shape of the ISW correla- 
tion, i.e. 



Pi.i = 



B(l)C\ 
; otherwise 



r ISW ; 2 < 2 < L 



(30) 



This definition implies that c\ directly measures the 
galaxy bias by optimally combining information from 

multipolcS < lmax- 

We are not restricted to estimating Cf T , but can si- 
multaneously measure Cf T and Cf 9 . Since doing so is 
not the goal of this paper, we simply measure the ampli- 
tude of the input Cf T and Cf 9 analogous to the template 
bandpower above. This provides a useful runtime check 
of our pipeline, as well as allowing us to estimate the 
galaxy bias. 

We estimate the c, by forming quadratic combinations 
of the data [I1G3, 



rx'c; 



'da 



diag 



(31) 



These are related to the estimated Cj by the Fisher ma- 
trix, F, 



- E( F-1 )* 



'a 



where 



Fa = -tr 

y 2 



■ ! dC ! dC 

diag 9c . diag 9c . 



(32) 



(33) 



If Cf -C y Cf g Cj T , then the Cj are a good approxima- 
tion to the maximum likelihood estimates of the a. The 
covariance matrix of the Cj is the inverse of the Fisher 
matrix, if the fiducial power spectra and noise used to 
compute Cj iag correctly describe the data. These as- 
sumptions must be tested and calibrated with simula- 
tions CSec. lVI A|) . 

Implementing the above algorithm is complicated by 
the sizes of the datasets; the data vector has 1,052,950 
elements making both storing and naively manipulating 
the covariance matrix impossible on presently available 
computers. Working in harmonic space is also not pos- 
sible due to the complicated geometry of the SDSS and 
WMAP masks. We implement the methods of |^ ex- 
tended to the sphere. All matrix-vector operations are 
performed using convolutions using the spherical har- 
monic transform code of 0, which scales as N 3 ^ 2 , com- 
pared to the iV 2 scaling of direct matrix multiplications. 
Matrix inversions are performed with a preconditioned 
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V. RESULTS 
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FIG. 4: The cross correlation between WMAP and the SDSS 
LRG density, as measured in flat bandpowers for the four 
WMAP bands considered in the paper. Also plotted is pre- 
dicted ISW signal for a !Jh = 0.3, Qa = 0.7 universe, scaled 
by the measured galaxy bias in Table|T] Note that these power 
spectra have been beam-deconvolved. 



Band 


bg {I < 50) bg {I < 400) 


Ka 


3.12 ± 1.65 4.20 ± 1.56 


Q 


2.77 ± 1.65 3.99 ±1.56 


V 


2.89 ± 1.65 4.05 ± 1.56 


w 


2.81 ± 1.66 3.91 ± 1.57 


V (flipped) 


0.39 ±1.58 


V (inverted) 


0.74 ± 1.57 



TABLE I: The measured galaxy bias as estimated from the 
ISW signal, for templates truncated at l ma x = 50 and l m ax = 
400. The errors on the data are computed from the Fisher 
matrix. Note that the different bands are correlated on the 
largest scales, explaining the similarities in the values of bias 
measured. Also shown are the measured biases for a flipped 
and inverted V band map; these maps are not expected to be 
correlated with the LRG density. 



The results of cross-correlating the WMAP tempera- 
ture maps with the the LRG density are shown in Fig.0] 
where the power spectrum has been estimated in flat 
bandpowers. For I < 50 where the ISW signal is ex- 
pected to be the strongest (Fig. 0) , we see evidence for a 
signal, peaking around I ~ 20. This signal is seen in all 
four bands and appears to be independent of frequency, 
consistent with the achromatic ISW effect. 

We want to estimate the amplitude of this signal, 
parametrized by the LRG bias. Performing a \ 2 fit to 
the above cross power spectra is undesirable for two rea- 
sons - (i) the flat bandpowers do not optimally use the 
information from the individual multipoles, and (ii) in- 
verting a noisy covariance matrix (due to the stochastic 
trace estimation) is likely to incorrectly estimate the sig- 
nificance of the detection. To avoid both these problems, 
we use the template bandpowers described in the pre- 
vious section to estimate the amplitude of the expected 
ISW signal (assuming our fiducial cosmology). We con- 



sider two values of L 



at which we truncate the ISW 



templates, l ma x = 50, and l max — 400. Truncating the 
template at large I ensures that we optimally use all the 
information present in the multipoles. However, we do 
see (Fig.QJ some large fluctuations at small scales, / > 50. 
As we discuss in the next section, these appear to be sta- 
tistical fluctuations, but choosing l max = 50 allows us to 
be conservative and ignore these multipoles. 

The results are in Table [I] and the best fit templates 
are plotted in Fig. 0] The signal again appears to be 
achromatic as expected. The values of the bias are con- 
sistent with measurements of the bias from the galaxy- 
galaxy power spectrum, b g = 1.82 ± 0.02, estimated as 
discussed in the previous section. We also note the val- 
ues of the bias obtained from truncating the templates 
at lmax = 50 and l m ax = 400 are consistent, although, 
not surprisingly, the error associated with l max — 50 is 
higher. 

Table [I] also presents the amplitude of the cross- 
correlation between the LRG density and flipped (north 
and south reversed) and inverted (antipodal points ex- 
changed) V band WMAP maps. These maps should not 
be correlated with the LRG density, and provide a useful 
systematic test internal to the data. Both biases obtained 
are consistent with zero as expected. 



VI. SYSTEM ATICS 



A. Pipeline Simulations 



conjugate gradient code 0; usm g the preconditioner 
in Appendix B of @; this typically converges to a frac- 
tional precision (with an L2 norm) of 10~ 8 in ~ 100 it- 
erations. Finally, the Fisher matrix is computed using 
a Z2 stochastic trace algorithm |4^; approximately 25 
random vectors achieves the necessary precision. 



End-to-end simulations are essential both to validate 
the pipeline, as well as to calibrate errors obtained from 
the Fisher matrix. Since realistically simulating the 
galaxy population would involve understanding the for- 
mation of LRGs and therefore, is not currently feasible, 
we simulate the microwave sky and cross-correlate these 
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Band 


b g (I < 50) bg (I < 400) 


X A (Sat) 


X A {Imax < 50) 


X 2 (Imax < 400) 








(16 dof) 


(8 dof) 


(2 dof) 


Ka 


0.30 ± 1.64 


0.30 ±1.56 


16.88 (0.39) 


9.90 (0.27) 


4.21 (0.12) 


Q 


0.30 ± 1.64 


0.29 ± 1.56 


19.29 (0.25) 


8.73 (0.37) 


5.78 (0.06) 


V 


0.29 ±1.65 


0.30 ±1.57 


15.00 (0.52) 


7.64 (0.47) 


3.99 (0.14) 


w 


0.28 ±1.69 


0.25 ± 1.60 


20.82 (0.19) 


8.51 (0.38) 


4.10 (0.13) 



TABLE II: The average galaxy bias derived from the 110 WMAP simulations cross-correlated with the actual LRG overdensity 
map, for l ma x = 50 and lmax = 400. Recall that the simulated microwave maps are uncorrelated with actual LRG map, but 
that the different WMAP frequencies are correlated. Also shown is the \ 2 °f the average power spectrum compared with zero 
and (in parentheses) the probability that \ 2 would be greater than the measured value for a correct model. Note that we also 
show the \ 2 value for the flat bandpowers in Fig. 2] 



realizations with the actual LRG density map. Since the 
ISW cross-correlation is much smaller than the individ- 
ual auto-power spectra, we simply generate microwave 
sky maps that are uncorrelated with the LRGs. 

The primary microwave temperature fluctuations are 
simulated by generating Gaussian random fields with 
an angular power spectrum identical to the best-fit 
ACDM spectrum determined by WMAP These 
maps are then convolved with the appropriate instru- 
ment beams [4^| : we assume circular beams for simplic- 
ity. Note that the temperature fluctuations are corre- 
lated between the different frequency bands, allowing us 
to estimate the frequency correlation of the power spec- 
trum estimates. We add noise by using the 110 simulated 
noise maps provided by the WMAP team [27j. In addi- 
tion to white noise, these maps also simulate the 1// 
detector noise and inter-detector correlations [51|. We 
make no attempt to add either Galactic or extragalactic 
foregrounds to these maps; these are small (especially for 
the V and W bands) outside of the KpO mask|5^|. All 
these maps are then masked with the Kp0nS10\ps mask, 
and then correlated with the actual LRG density in an 
identical manner to the actual temperature data. 

The results from the 110 simulations are summarized 
in Table[n] The frequency correlations on large scales are 
evident from the identical values of the average bias ob- 
tained. Note that measured power spectra are consistent 
with zero, both for the flat and template bandpowers, 
with x 2 /dof ~ 1 in all cases. Finally, we observe that 
the error on the bias from the Fisher matrix is consistent 
with the run-to-run simulation error (within the variance 
of the simulations, Aa/a ~ l/2\/55). 



B. Galactic Foregrounds 

The two sources of contamination in the galaxy cata- 
log that could correlate with Galactic microwave emission 
are stellar contamination and incorrect Galactic extinc- 
tion corrections. The Kp0nS10\ps mask by construction 
excludes regions with the worst contamination; we test 
for any residual contamination by directly estimating the 
level of contamination (Sec. lVI B 1(1 and cross-correlating 
foreground emission templates with the galaxy density 
(Sec. lVlB2t . 



1. Stellar density and Galactic extinction contamination 

In order to estimate the level of stellar contamina- 
tion and incorrect extinction corrections, we compute the 
zero-lag correlation between the LRG density maps and 
maps of the stellar density and Galactic extinction. To 
avoid discretization effects, we re-pixelize the LRG cat- 
alog onto a HEALPIX sphere with 49152 (resolution 6) 
pixels. We estimate the stellar density by selecting stars 
identified by the SDSS photometric pipeline, with r band 
PSF magnitudes between 18.0 and 19.5, and pixelizing 
these onto the same HEALPIX sphere. The extinction 
map is constructed directly from the E(B — V) reddening 
map of |53(. 

Contour plots of the LRG density as a function of 
stellar density and Galactic extinction are in Fig. |S| 
the contours show conditional probabilities of 5%, 50% 
and 95%. The near-horizontal contours (even at the ex- 
tremes) suggest negligible levels of contamination. The 
cross-correlation coefficients of r = 0.006 (stellar density) 
and r = 0.004 (extinction) further support this conclu- 
sion. 



2. Foreground Templates 

A second test of Galactic foreground contamination is 
to cross-correlate templates of Galactic emission with the 
LRG density maps. We construct these templates with 
three components: 

• Thermal dust emission : We use the 100 micron 
maps of dust emission [EiJ, extrapolated to the 
WMAP frequencies using the two-component dust 
model (Model 8) of H3- 

• Free-Free emission : We model the free-free com- 
ponent with the all sky Ha maps of [55j using the 
prescription in [52] . 

• Spinning Dust/Synchrotron Emission: 
WMAP shows evidence for an additional source 
of emission at low frequencies due either to syn- 
chrotron emission |52| , or spinning dust |56j . Since 
the source of this emission is still under debate, we 
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FIG. 5: Conditional contour plots of the galaxy density as a 
function of the stellar density and galactic extinction, E(B — 
V). The contours are the 5%, 50% and 95% contours, while 
the horizontal line is the mean galaxy density. Also shown 
are the correlation coefficients, r for both datasets. There is 
no evidence for significant stellar contamination or incorrect 
extinction corrections in these data. 



FIG. 7: The upper panel shows the redshift distribution in- 
ferred from two different levels of regularization. The solid 
line is the distribution used throughout this paper, while the 
dashed line is less regularized by a factor of 10. The lower 
panel compares the templates for Qm = 0.2, 0.3, 0.4 [short 
dashed, solid, long dashed] universes using the above redshift 
distributions. 
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adopt the simple phenomenological prescription in 
56] to model this component. 

These maps are constructed in the same pixelization as 
the WMAP temperature maps, and are masked with the 
Kp0nS10\ps mask. The resulting maps (identical in form 
to the temperature maps) are then analyzed to compute 
the cross-correlation power spectrum. 

The results are shown in Fig. EI As expected, the con- 
tamination increases with decreasing frequency, with Ka 
showing the worst contamination and V and W with neg- 
ligible contamination. However, Galactic foreground con- 
tamination is no more than 20% of the detected signal 
for Ka, and is less than 10% in the other bands. There- 
fore, while Galactic foregrounds may be responsible for 
the slightly higher value of the bias in the Ka band (Ta- 
ble [U, they cannot explain most of the signal in the 
temperature-galaxy correlation. Unlike the ISW effect, 
extragalactic foregrounds, particularly point sources, are 
expected to show a dependence on frequency. Tableland 
Fig.0]suggest that the any such frequency dependence is 
much less than the errors on the measurements. 



Redshift distribution uncertainties 



FIG. 6: The cross correlation between the WMAP foreground 
templates and the SDSS LRG density, measured in the same 
bandpowers as Fig. 0] The contamination from known fore- 
grounds is clearly subdominant to the CMB temperature cross 
correlation. 



Predicting the ISW signal requires knowing the red- 
shift distribution of the LRGs. As discussed in 0], es- 
timating the redshift distribution from measured photo- 
metric redshifts involves deconvolving the effect of photo- 
metric redshift errors. This process is unstable and must 



10 



* 1 

OJ 

o 

P 
+ 



10 



100 



1000 



FIG. 8: The cross correlation between the WMAP V band 
and the LRG density, measured for the two contiguous SDSS 
regions separately. The measured power spectra are consis- 
tent within the estimated errors. Note that the bands are ~ 
twice the size of bands in Fig. 2] and the dashed points have 
been displaced for clarity. 



be regularized; unfortunately the choice of the regular- 
ization parameter and therefore, the estimated rcdshift 
distribution, is not unique. Fig. [7| estimates the error 
that this uncertainty introduces into the measurement of 
the bias. The upper panel compares the redshift distri- 
bution used throughout this paper with a distribution 
estimated with a smaller regularization parameter. The 
principal effect of reducing the regularization is ringing 
in the inversion, caused by noise in the photometric red- 
shift distribution. However, as is evident from the lower 
panel, changing the redshift distribution makes only sub- 
percent changes to the ISW signal. This insensitivity de- 
rives from the fact that the greatest contribution to the 
ISW signal is coming from z ~ 0.5, where the inversion 
is most stable. In addition, integrating along the line of 
sight makes the signal principally sensitive to broadband 
features in the redshift distribution. 



D. Excess power at 



50, 100, 200 



The measured cross correlation (Fig. shows a 
marginally significant (~ 2a) excess at multipoles around 
Z — 50, Z — 100 and / ~ 200. This could either be a sta- 
tistical fluctuation, or it could represent foreground con- 
tamination. As discussed above, there are a number of 
possible systematics - we consider each of these in turn 
below. 
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FIG. 9: The cross correlation between the WMAP V band 
and the LRG density, measured in Az = 0.1 photometric red- 
shift slices. The bands are exactly those used in Fig. 0] Note 
that the power spectra are consistent within the estimated er- 
rors. Also plotted is the ISW template scaled by the measured 
V band bias. 



(Fig.jBJ, this would either require an order of mag- 
nitude mis-estimate of the known foregrounds, or 
a new component. Combined with the fact that 
this would have to be significant at high Galac- 
tic latitudes only (since the known templates agree 
well with the measured microwave emission in all 
regions outside the KpO mask), Galactic contami- 
nation is unlikely to be the source of this power. 

• Extragalactic foregrounds : Extragalactic fore- 
grounds broadly divide into two classes - SZ con- 
tamination and microwave point sources. Known 
point sources have very strong frequency dependent 
spectra, not observed in Fig. 0] The SZ effect is 
characterized by a temperature decrement |57) at 
these frequencies, the opposite of what is observed. 
The kinetic SZ effect, although frequency indepen- 
dent, is significantly smaller than the thermal SZ 
effect and should be negligible. 

• Underestimation of Errors : A third possibility is 
that the errors (derived from the Fisher matrix) are 
underestimated. However, comparing with simula- 
tions, we find that the Fisher errors agree to better 
than - 10%. 



• Galactic contamination : Given the size of this ef- Since known systematics are unlikely to cause this ex- 

fect and the expected contamination on these scales cess power, we divide the LRG catalog into various sub- 
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FIG. 10: The estimated galaxy bias from the galaxy autocor- 
relation compared with the bias (and error) estimated from 
the ISW signal (shaded region), as a function of Qm (assum- 
ing Qm + Qa = 1). The errors on the autocorrelation bias are 
negligible compared to the ISW errors. The best fit value of 
Q M is Q M = 0.20±JJ;J5|(lff)±g;S§(2<7). Note that the likelihood 
for Qm is extremely non- Gaussian (see Table ITTO . 



samples to attempt to localize this effect. Fig. [S] shows 
the cross-correlation with the two contiguous regions of 
the SDSS (see Fig- EJ, while Fig. cross-correlates dif- 
ferent photometric redshift slices. Within the errors, 
these different subsamples are consistent with each other 
around I ~ 50; the redshift bin from z = 0.3 to 0.4 shows 
an excess of power at I ~ 100. While this excess might be 
caused by a failure of the LRG selection algorithm at the 
junction of Cut I and Cut II, a precise mechanism that 
would cause an excess in power when correlating with the 
CMB is not apparent. 

At the current level of accuracy, it is impossible to 
determine if this excess is simply statistical, or if it in- 
dicates systematic effects. Greater sky coverage as the 
SDSS nears completion could help answer this question. 
We note that this excess signal adds to the best fitted 
value of ISW amplitude in our analysis (compare columns 
2 and 3 of Table UJ) , but at a small level since these small 
scale modes are downweighted in the theoretical tem- 
plates (Fig. (IJ. This differs from other analyses where 
this procedure was not adopted 01 ■ I n addition, most 
of previous work has been based on correlation function 
analyses, where correlations between the bins in the cor- 
relation function make a scale dependence of the signal 
difficult to identify. 



VII. COSMOLOGICAL CONSTRAINTS 

The fact that the ISW signal — > oo as Qm — > suggests 
a method to place a lower bound on the matter density 
within the context of models with Qm + = 1. If one 
computes the bias both from the galaxy autocorrelation 
function (b g ) and the ISW effect (bisw), then one finds 
that bisw — > as Qm — > faster than b g , as does the 
error on bisw- Comparing the two values of the bias 
allows one to constrain the value of fi m- Note that this 
only puts a strong lower bound on Qm- The upper bound 
is determined by how strongly one detects the ISW effect; 
Qm = 1 implies no ISW effect. Therefore, one would 
only expect to be able to rule out Qm — 1 at ~ 2.5er, our 
quoted detection significance. 

The results of performing this exercise are summarized 
in Fig, mil and Table ITTT1 The galaxy and ISW bias are es- 
timated by using their predicted templates from Sec. [H] 
We have limited ourselves to the WMAP V band both 
for computational convenience and because the effects of 
galactic foregrounds are minimum there. We continue to 
use the fits to the measured power spectra used in previ- 
ous sections to compute C~ l as these correctly measure 
the cosmic variance of these data. 

The best fit value for the matter density is Qm = 0.20 
with lcr limits of Qm = 0.17 and Qm = 0.26. The 
strength of the lower bound is evident from Table IIIII 
where we see that Qm = 0.15 is ruled out at ~ 2er, while 
Qm = 0.10 is ruled out at > 8er. Modest increases in 
the ISW detection significance will translate into stronger 
constraints. For instance, increasing the survey area by a 
factor of 2 (e.g., with the SDSS area on completion) could 
rule out Qm = 0.15 at > 3er while an all sky survey at 
moderate redshift could rule it out at > 6o\ While these 
bounds will not replace traditional parameter estimates 
(for eg. 58] ) , they will provide much needed independent 
tests of dark energy. 

So far we have held all other cosmological parameters 
fixed in the above analysis. Marginalizing over these pa- 
rameters, in general, will weaken these constraints; these 
however will only weakly change the strength of the lower 
bound. For example, reducing og by a factor of x in- 
creases the ISW bias by a factor of x 2 , but the galaxy- 
galaxy bias only by a factor of x. However, decreasing erg 
from 0.9 to 0.8 (a 10% reduction) reduces the significance 
of Qm = 0.15 to 1.6cr from 2.0c, while making erg = 1 
increases it to 2.5<r. 

The above discussion assumes that the acceleration is 
due to a cosmological constant, i.e. due to a component 
with equation of state, P = —p. An important general- 
ization are quintessence models with equation of state, 
P = wp. For w > — 1, the ISW signal is enhanced, in- 
creasing the best fit value of Qm relative to the w = — 1 
case, while w < — 1 reduces the signal, allowing for even 
lower values of Qm', the likelihoods for Qm assuming dif- 
ferent values of w are shown in Fig. 1111 Unfortunately, 
the weak ISW detection does not, by itself, allow us to 
place useful constraints on Qm and w jointly. However, 
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1.31 


0.24 


0.13 


8.3 


7.2 


9.5 


0.26 


0.13 


8.2 


7.1 


9.4 


0.15 


1.62 


0.74 


0.38 


2.3 


1.8 


2.8 


0.85 


0.38 


2.0 


1.6 


2.5 


0.20 


1.73 


1.38 


0.73 


0.5 


0.2 


0.8 


1.71 


0.71 


0.0 


0.2 


0.3 


0.25 


1.79 


2.09 


1.15 


0.3 


0.4 


0.1 


2.79 


1.11 


0.9 


1.1 


0.7 


0.30 


1.82 


2.89 


1.66 


0.6 


0.8 


0.5 


4.06 


1.57 


1.4 


1.6 


1.3 


0.35 


1.83 


3.79 


2.27 


0.9 


1.0 


0.8 


5.55 


2.09 


1.8 


1.9 


1.7 


0.40 


1.85 


4.81 


2.98 


1.0 


1.1 


0.9 


7.27 


2.69 


2.0 


2.1 


1.9 


0.45 


1.85 


5.99 


3.84 


1.1 


1.1 


1.0 


9.28 


3.37 


2.2 


2.3 


2.1 


0.50 


1.85 


7.39 


4.88 


1.1 


1.2 


1.1 


11.65 


4.18 


2.3 


2.4 


2.3 


0.55 


1.85 


9.04 


6.14 


1.2 


1.2 


1.1 


14.49 


5.14 


2.5 


2.5 


2.4 


0.60 


1.85 


11.09 


7.72 


1.2 


1.2 


1.2 


17.97 


6.32 


2.6 


2.6 


2.5 



TABLE III: The estimated bias from the galaxy autocorrelation (b g ) and the ISW cross-correlation (bisw) as a function of 
Qm, shown for both templates truncated at l max = 50 and l m ax = 400. The significance, [a] , is defined as \b g — bisw\/o~b,isw 
where &b,isw is the error on the bias inferred from the ISW effect, and Ob,g 1- This describes number of a a particular 
value of Qm is away from the best fit value and can be used in likelihood analyses of cosmological parameters (\ 2 = [°] 2 )- We 
also show the effect of changing as on these results; increasing as increases the inferred value of Qm- The ISW effect is seen 
to provide a strong lower bound on the value of Qm. 
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FIG. 11: The likelihood for Qm for different values of the 
equation of state, w, computed using the prescription of Ap- 
pendix^ The solid line assumes w — —1, i.e. a cosmological 
constant, while the dashed lines assume (from left to right) 
w = —1.2, —0.8, and —0.5 respectively. The crosses show 
the likelihood computed using the full quadratic estimator 
discussed above, for w = —0.8, and test the approximations 
made in the appendix. The dotted lines denote the 1-, 2-, and 
3-a levels. 



one can combine the ISW constraints with other data. 
This would constrain models with w < — 1 that require a 
low Qm, currently disfavored by other data. We present 
an simple algorithm for computing the likelihood for cos- 
mological models parametrized by (Qm,o~s,w) from the 
ISW signal in Appendix lAl 

Finally, we make a parenthetical observation about the 
size of the errors for templates truncated at different l m ax 



values. We observe that for low Qm universes, truncating 
at different multipoles has little effect, while for higher 
f)jfi there is a significant difference. This is a direct re- 
sult of the template weighting; at low Qm, most of the 
weight is at low multipoles where the signal peaks and so, 
including higher multipoles has little effect. As Qm in- 
creases, the higher multipoles gain importance, decreas- 
ing the error. 



VIII. DISCUSSION 

We have cross-correlated the microwave temperature 
maps observed by WMAP with LRG overdensities ob- 
tained from the SDSS. Our goal is to extract the maxi- 
mum amount of information on ISW available from the 
data. The major differences between our analysis and 
previous ISW analyses are: 

• Power Spectrum: Most analyses (with the excep- 
tion of [fij) have used the correlation function in- 
stead of the power spectrum. While it is true that 
theoretically the power spectrum and correlation 
function are just Legendre transforms of each other, 
this is no longer true in the presence of noise and 
sky cuts. In particular, the correlation function has 
extremely correlated errors, while the power spec- 
trum errors are almost uncorrelated on large scales. 
The presence of correlated errors makes an accurate 
determination of the covariance matrix with either 
Monte Carlo or jack-knife techniques more difficult. 

• Optimal C _1 Weighting: We weight our data 
with the inverse covariance, resulting in minimum- 
variance error bars even in the presence of sky cuts. 

• Fisher Matrix errors: We estimate our errors using 
the Fisher information matrix, instead of jack-knife 
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error estimates favoured by a number of other anal- 
yses. Jack-knife errors tend to underestimate ISW 
errors because of the small number of uncorrelated 
jack-knife patches on the sky. These small numbers 
of samples also results in noisy covariance matrices, 
potentially biasing \ 2 fitting. We also note that 
the Fisher matrix errors are the smallest possible 
errors, given the intrinsic sample variance. We ver- 
ify that we saturate this bound with simulations, 
indicating an optimal measurement. 

• x 2 Fitting: All analyses previous to this one have 
estimated the significance of their detection by fit- 
ting power spectrum or correlation function bins 
to ISW templates. We avoid this step by directly 
estimating the amplitude of a fiducial template; 
the choice of this fiducial template affects our con- 
clusions about the significance of the effect only 
marginally, as seen in Table ITTT1 This has two ad- 
vantages - (i) it optimally uses the information from 
all multipoles, and (ii) it avoids systematics intro- 
duced by x 2 fitting with a noisy covariance ma- 
trix. The noise in the covariance matrix is further 
exacerbated by strong correlations between bins; 
this results in badly conditioned covariance matri- 
ces where errors in the most singular directions (the 
ones with the smallest eigenvalues) will dominate 
the fits. Note that these problems are made worse 
with the use of jack-knife errors, due to correlations 
in the different jack-knife samples on large scales 
and the small number of samples available. 

The two analyses that use galaxies samples closest to 
ours are and 0]; we n0 ^ e that both these analyses 
yield results similar to ours within the quoted errors. 
However, in addition to the discussion above, our analy- 
ses differ in a number of areas. The most trivial is in the 
difference in the sky coverage, ~ 4000 deg 2 compared 
to the ~ 2000 and ~ 3400 deg 2 respectively. Further- 
more, we use accurately calibrated photometric redshift 
distributions, including a deconvolution of the redshift er- 
rors 0[ . Most importantly, we use theoretical templates 
to optimally weight the different scales, immunizing us 
to contamination from statistical fluctuations on smaller 
scales. 

Despite our attempt to be nearly optimal and our use 
of the latest and largest SDSS data sample we do not 
find a strong detection of ISW. We observe a correlation 
on large scales at 2.5 a, consistent with the ISW effect. 
We fit to the amplitude of this effect, optimally combin- 
ing the information from multipoles < 400 by using a 
predicted template, and, for CIm = 0.3 and as = 0.9, 
obtain a bias for the LRGs of b g = 4.05 ± 1.56 for the V 
band correlation, and similar values for the other bands. 
Restricting to multipoles < 50 yields b g = 2.89 ± 1.65, 
consistent with the previous value although with a some- 
what larger error. These are consistent with the value 
from the galaxy-galaxy auto-correlation, b g = 1.82 ±0.02 
and even better agreement is obtained for a smaller value 



of n M (~ 0.2). 

We explore systematic effects that could contami- 
nate the signal. Known Galactic and extragalactic fore- 
grounds are subdominant to the measured ISW signal. 
However, we do detect a marginally significant (~ 2a) 
correlation at scales of I ~ 50 and I ~ 100. The frequency 
dependence and sign of the correlation are inconsistent 
with extragalactic foregrounds, and Galactic foregrounds 
are approximately an order of magnitude too small to in- 
fluence our results. Subdividing the LRG sample both 
on the sky, and in redshift slices, yields power spectra 
consistent with each other, further ruling out systematic 
effects. We tentatively conclude that this excess power 
is a statistical fluctuation, and wait for future data to 
confirm or disprove this conclusion. This small scale ex- 
cess power slightly increases the statistical significance of 
ISW in our analysis, but less so than in previous analyses 
where theoretical template weighting was not used. 

The ISW effect provides a useful and independent 
probe of the properties of the dark energy. Although the 
weakness of the current detections preclude detailed pa- 
rameter estimations, we attempt to estimate the matter 
density, by comparing the ISW bias to the value obtained 
from the auto-power spectrum analysis. The rapidly in- 
creasing amplitude of the ISW signal with decreasing 
matter density provides a strong lower bound on the mat- 
ter density, ruling out VLm = 0.15 at 2a. These represent 
fits to Q,m only and assume ct 8 = 0.9 and w = — l|5S|. 
The current ISW detection is not significant enough to 
jointly constrain CIm and w. However, these constraints 
can be incorporated into parameter estimation efforts; we 
provide a simple likelihood prescription to do this. As 
an example, assuming w = —1.3 reduces the predicted 
ISW signal and the matter density must be decreased to 
match observations: we find the best fit value for matter 
density is reduced to f2 m = 0.14 with a 2a upper limit 
of 0.29. Such a low value is in conflict with other deter- 
minations of the matter density (see e.g. [59)) and as a 
result w < — 1 models are strongly constrained. 

In the present analysis we have ignored any evolution 
of the bias with redshift, since the redshift distribution 
is relatively narrow. Furthermore, estimates of the bias 
from the galaxy-galaxy power spectrum for different red- 
shift slices yields a bias of between 1.7 and 1.9; this vari- 
ation is too small to affect our analysis given the current 
errors. However, this could be a more significant issue 
for heterogeneous catalogs like the X-ray background and 
the NVSS radio catalog. 

ISW detection efforts are still in their infancy, and 
future observations should lead to a stronger detection. 
This will hopefully realize the full potential of using the 
ISW effect as a probe of dark energy. 
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0.41 
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1.99 
2.13 



0.15 
0.27 
0.41 
0.58 
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1.00 
1.25 
1.51 
1.77 
2.00 
2.16 



0.15 
0.27 
0.42 
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1.00 
1.24 
1.50 
1.76 
1.99 
2.16 



0.16 
0.28 
0.42 
0.59 
0.78 
1.00 
1.24 
1.49 
1.74 
1.97 
2.15 



TABLE IV: Approximate multiplicative factors that scale 
the ISW power spectrum as a function of w and Q,m. The 
bias scales as biswiw) = bisw(w = — l)/wti as . Note that 
the fiju dependence is weak, and can be ignored at the level 
of accuracy required. 
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where the right hand side asssumes as = 0.9, and w — 
— 1, i.e. the values in Table ITTT1 



The effect of changing w on b g can be approximated 
by rescaling with the ratio of the growth factor D(z), 



APPENDIX A: COSMOLOGICAL PARAMETER 
ESTIMATION 

We discuss a prescription to extend the results of Ta- 
ble [ffi] to compute the likelihood of a cosmological model 
parametrized by (Om, os, w). The likelihood is computed 
with a x 2 function, comparing the galaxy bias obtained 
from the autocorrelation and ISW analyses, 



X" 



b g (n M ,<78,w) - b IS w(^M,(J8,w) 



<^b,ISw(^M, 0"8, W ) 



(Al) 



As in Sec. IVIll we ignore the error on b g , since it is sub- 
dominant to the error on the ISW measurement. 

Table IIIII presents measurements of b g , bisw, and 
Cbjsw as a function of 17m, for as = 0.9, and w = —1. 
Extending the results for different as is straightforward, 
as discussed in Sec. I VIII 
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where z mean ~ 0.5 is the median redshift of the LRG 
catalog. Computing the effect on bisw would, in gen- 
eral, require a recomputation of the predicted ISW sig- 
nal; however, as Fig. 1121 indicates, we can approximate 
this (at the level of accuracy demanded by the strength 
of the ISW signal) by a simple multiplicative bias Wbias , 
tabulated in Table llVl The dependence of Wbias is weak, 
and can be ignored. Therefore, we can complete our spec- 
ification of bisw and a^jsw with 



bisw(^M,as,w) 

&b,ISw(QM,&8,w) = 
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FIG. 12: The ISW angular power spectrum as a function of w, 
starting at the bottom with w = —1.3 and increasing in steps 
of 0.1, relative to w = — 1. The cosmological parameters used 
are Q M = 0.3, fl A = 0.7, fl b = 0.05 and h = 0.7. The power 
spectra can be well approximated by a simple multiplicative 
bias. 
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